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ABSTRACT 

Scale transformations have played an extremely successful role in studies of cosmolog- 
ical large-scale structure by relating the non-linear spectrum of cosmological density 
fluctuations to the linear primordial power at longer wavelengths. Here we general- 
ize this approach to investigate the usefulness of scale transformations for nonlinear 
higher-order statistics, specifically the bispectrum. We find that the bispectrum pre- 
dicted by perturbation theory at tree-level can be rescaled to match the results of full 
numerical simulations in the weakly and intermediately nonlinear regimes, especially 
at high redshifts, with an accuracy that is surprising given the simplicity of the pro- 
cedure used. This discovery not only offers a simple practical way of calculating the 
matter bispectrum, but also suggests that scale transformations may yet yield even 
deeper insights into the physics of hierarchical clustering. 
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1 INTRODUCTION 

In the linear stage of gravitational instability, cosmological 
density perturbations evolve in such a way that each Fourier 
mode grows at the same rate as the others, preserving the 
shape of the primordial power spectrum. In the later stages, 
things become much more complicated, with mode-mode 
interactions changing the relative growth rates of different 
harmonics and coupling their phases together to generate 
non-zero higher-order spectra (Coles & Chiang 2000). Nev- 
ertheless, inspired by the creative work of Hamilton et al. 
(1991), empirical formulae have been established that can 
accurately predict the two-point correlation function ^(r) or 
power spectrum P{k) of cosmological density fluctuations 
in the non-linear regime given only the initial linear power 
spectrum and the background cosmological model (Jain, Mo 
& White 1995; Peacock & Dodds 1996). The key element of 
these formulae is a scale transformation which expresses the 
conservation of particle pairs which, in turn, is one of the in- 
finite set of equations that forms the BBGKY hierarchy, c.f. 
Peebles (1980). The discovery of this unexpectedly success- 
ful transformation for the nonlinear ^(r) or P{k) was a great 
leap forward in the understanding of cosmological clustering 
evolution, but with the renaissance of halo model (Cooray 
& Sheth 2002), and the emergence of various high-level per- 
turbative techniques (e.g. Szapudi & Kaiser 2003; McDonald 
2007), interest in the classic form introduced by Hamilton et 
al. (1991) has to some extent faded away. Nowadays it seems 



much more fashionable to calibrate gravitational nonlinear- 
ity by tuning parameters of the halo model or using empir- 
ical fitting formulas based on it (Smith et al. 2003). 

After great successes with two-point correlation func- 
tions, it was widely anticipated that the halo model would 
be able to generate precise nonlinear three-point statistics 
over a full range of length scales (Scoccimarro et al. 2001; 
Takada & Jain 2003). Detailed examination against simu- 
lations, however, has proved somewhat disappointing. New 
elements have had to be incorporated to improve the current 
halo model even after introducing ad hoc free parameters ac- 
counting for mass cutoff and halo boundary (Fosalba, Pan 
& Szapudi 2005), or replacing spherical haloes with triaxial 
ones (Smith, Watts & Sheth 2006). 

In this paper, we revive the scale transformation idea 
to derive a model that approximates the nonlinear contri- 
butions to the bispectrum of cosmological matter pertur- 
bations. The bispectrum is the lowest order (and therefore 
the simplest) diagnostic of non-Gaussianity because it van- 
ishes identically for any Gaussian random field regardless 
of its power spectrum. If the initial density fluctuations are 
Gaussian, as is expected in the simplest inflationary models, 
the bispectrum is therefore completely induced by non-linear 
gravitational evolution. However, these higher-order effects 
on the bispectrum are difficult to calculate with reasonable 
accuracy even on large scales where the evolution is quasi- 
linear. Differences between full numerical simulations and 
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the predictions of Eulerian perturbation theory at tree-level 
exceed ~ 10 — 20% at a scale k ~ O.l/i/ Mpc. The pertur- 
bation theory at one-loop level, the next possible improve- 
ment on the tree-level theory, requires heavy numerical in- 
tegration and has limited success (Scoccimarro et al. 1998). 
Lagrangian perturbation theory at second- and third-order 
provides a better template than Eulerian theory, but its dy- 
namic range is heavily restricted by the shell-crossing scale 
at around k ~ OAh/ Mpc (Scoccimarro 2000). A further in- 
convenience of adopting Lagrangian perturbation theory is 
that the bispectrum it predicts is not computed, but mea- 
sured from N-body simulations of which the particle motions 
are controlled by Lagrangian theory. 

One practical way forward is to modify the kernel of 
perturbation theory at tree-level to develop empirical fit- 
ting formulae such as extended and the hyper-extended per- 
turbation theory (Colombi et al. 1997; Scoccimarro & Frie- 
man 1999). However, the performance of these techniques 
known to be quite poor, even in cases where the initial 
power-spectrum is completely scale-free (Hou et al. 2005). 
The best results so far in this vein relate to the empirical 
model of Scoccimarro & Couchman (2001, herefter SC2001), 
using hyper-extended perturbation theory, but its average 
deviation from A'^-body results is at the level of 15%, still 
way beyond the accuracy levels expected in the era of pre- 
cision cosmology. On the other hand, even this produces 
some dubious features, such as the apparently trough in the 
bispectrum at scales ~ 0.03 < k <~ 0.12/i/ Mpc, where 
baryon acoustic oscillations leave their footprint on large- 
scale structure statistics. 

An accurate template for the bispectrum on large scales 
is extremely important in the epoch where large galaxy red- 
shift surveys are available to unveil details of large scale 
structure and to constraint cosmological parameters with 
higher-order statistics. It is possible that an accurate model 
can be achieved based on the formula of SC2001 facilitated 
with high precision simulations in huge box, or perhaps with 
a theory established by the renormalization technique of Mc- 
Donald (2007). However, in this paper we will show that, at 
least in the quasi- and intermediate nonlinear regime, it is 
possible to recover the bispectrum quite accurately using the 
scale transformation described above, an approach which is 
free from the laborious calibration of fitting parameters and 
formidable multi-dimensional integrations. 



2 SCALE TRANSFORMATIONS AND THE 
BISPECTRUM 

2.1 Scale Transformations of Second-order 
Statistics 

Let ^(r) be the two point correlation function in real space at 
scale r. The number of neighbours of a central point within a 
spherical top-hat window of radius r is given by '"^[1 + C('")]5 
where ^(r) — 3 Jq s'^£,{s)ds/r'^ is the volume- averaged two 
point correlation function. Hamilton et al. (1991) proposed 
that by pair conservation at any stage of gravitational evo- 
lution, one can quote a Lagrangian separation vl defined 
by 

ri =r^L [I+Htnl)] , (1) 



at which the nonlinear function ^j^j^ is an universal function 
of the linear one: ^j\il{i~nl) ~ /[CL(''i)l- The effectiveness of 
the paradigm is at least partly explained by Nityananda & 
Padmanabhan (1994): such a nonlinear scaling relation can 
arise if the pairwise velocity demonstrates certain scaling 
features. In fact, based on the universal properties of the 
pairwise velocity distribution, the pair conservation equation 
can be solved with an iterative technique to produce the 
correct nonlinear ^ from the input linear function (Caldwell 
et al. 2001). 

In Fourier space, the same sort of scaling can be applied 
to the dimensionless power spectrum A^(fc) = P{k)k'^ / 
by a relation of the form 

fcz. = [1 + ^l^ikML)] kNL , (2) 

so that the nonlinear power spectrum PNL(kNL) is obtained 
through A^^(fcjvi) — f [A|^(fc_L)] . The scale transformation 
of Eq. ([2| does not have such a firm theoretical motivation 
as its real space version, although there is no question that 
works very well in practice (Peacock & Dodds 1996). 

2.2 The Bispectrum at Tree- level 

The bispectrum is the three point correlation function de- 
fined in Fourier space. Denote the Fourier transformation of 
the cosmic density contrast with 5(k). The bispectrum is 

B(ki,k2,k3) = (5(ki)5(k2)r(k3)fc(ki+k2+k3)) , (3) 

in which Sn is the Dirac function. In a statistically 
isotropic universe, i3(ki,k2,k3 = — ki — k2) is written as 
_B(fci, fc2, fcs), with ka = \/kf + k2 + 2kik2fJ,i2 and /ii2 = 
ki ■ k2/(feiA;2). 

Because the initial distribution of dark matter is Gaus- 
sian, there is no non-zero bispectrum in the scheme of linear 
evolution. The first non-trivial contribution to bispectrum, 
at tree-level, comes from the second order terms in the ex- 
pansion of the cosmic density fiuctuation. Explicitly, 

BpT{ki,k2,ka)=PLiki)PL{k2)F2{ki,k2,fii2)+cyc. , (4) 

where eye. refers to the other two terms obtained by making 
cyclic permutations of the indices of the first term, Pl is the 
linear power spectrum, and F2 is the kernel of second order 
in Eulerian perturbation theory (Goroff et al. 1986): 

f2(fci,fc2,Mi2) = y + M12 + + ^M?2 • (5) 

Note the dependence on cosmological parameters of F2 is 
extremely weak, so it will be ignored here. 

2.3 Configuration Re-arrangement 

The three wave vectors defining the bispectrum form a trian- 
gular configuration, and it is well known that there is strong 
dependence of bispectrum on the shape of this triangle. If 
we assume that the transformation rule defined by Eq. ((2)1 
holds for arbitrary k, so that a new set of scales k in Fourier 
space is given by 

fcl,2,3 = [1 + A^i(fci,2,3)] fci,2,3 , (6) 

with A^^ computed by the formula of Smith et al. (2003), 
the shape of the original triangle is obviously not preserved. 
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Figure 1. Bispectrum of equilateral triangles, PT means the second-order perturbation theory, SC2001 is the prediction of Scoccimarro 
& Couchman (2001), PT rescaled is the expectation given by Eq. (10). The spurious trough at large scales of SC2001 can be seen in the 
z=0 panel. 



By the requirement of ki,2,3 ~ 0, the transformed triangle 
has new set of cosines Jl, e.g. 

M12 = ~(kl + fci - kl)/{2kik2) ■ (7) 

Now comes the crucial step in our model. We assume that 
the non-linear bispectrum is directly related to the tree-level 
bispectrum of the new configuration in Lagrangian space. 
In the non-linear regime this relationship is expected to be 
expressed by a very complicated functional form, but here 
we simply conjecture that we can use the linear expression, 

B{kx,k2, H12) oc BpT{ki,k2,TL\2) ■ (8) 

The remaining ingredient we need is the factor by which the 
amplitude of the bispectrum is boosted during the rescaling. 

2.4 Amplitude Boost 

We have tested many schemes to calculate the boost to be 
applied to the amplitude, but the best is one involves the 
idea that a non-linear fluctuation in Lagrangian space is 



expressible as a convolution of the linear density contrast 
with some filter so that, in Fourier space, 

&NL{kL) = w{kL)SL{kL) , (9) 

where ixj(fc) is the window function (the Fourier transform 
of the filter) and fci is given by Eq. ((2}. The form of this 
convolution is motivated by terms that arise in halo model 
from cross-correlation terms in the power spectrum and bis- 
pectrum. In this scenario, the nonlinear growth of the fiuc- 
tuation in a region is imagined to consist of two stages: (1) 
aggregation of clustering power in Lagrangian space from 
other regions via convolution; and (2) translation from the 
Lagrangian (initial) scale to the Eulerian (final) scale. 

Putting all the elements of Eq. (4) - (9) together, the 
expected nonlinear bispectrum in the recipe is 

_B(fci, fe2, M12) = w{ki)w{k2)w{k^)BpT{ki,k2,Jii2) 

r ~ ~ Tl/2 (10) 

= \PNL{h)/PL{h)^ , i= 1,2,3 . 

We have to address the issue that, in the model, the pos- 
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Figure 2. B{ki,k2,0) at z=0, (9 = cos ^ /ii2, see Fig. (1) for legend and labels. 



sible phase shifts induced by convolution are not taken 
into account. This might have a significant effect in the 
strongly non-linear regime, since the bispectrum is a phase- 
dependent function (Watts & Coles 2003). 



3 COMPARISON WITH SIMULATIONS 

Four outputs at 2 = 0, 0.5, 1, 2 of the VLS (Very Large Sim- 
ulation) data provided by the Virgo consortium are used to 
check our model. The simulation runs consist of 512'' parti- 
cles in a cubic box of size 479h~^ Mpc, and the cosmologi- 
cal parameters adopted are Qm = 0.3, flv = 0.7, F — 0.21, 
h = 0.7 and erg = 0.9 (MacFarland et al. 1998). In order 
to estimate cosmic variance, the VLS simulation is divided 
into eight distinct sub-cubes of half of the original size. The 
bispectrum is measured using the method of Scoccimarro et 
al (1998), and the Icr dispersion among the eight measure- 
ments is taken as an estimate of the error bars. 

To avoid any bias involved with taking ratios of two 
statistical quantities and leakage of the power spectrum into 
the bispectrum, we work with the bispectrum itself rather 



than the reduced bispectrum Qs 
throughout this paper. 



B/(PlP2+PlP3+P2P3 



3.1 Equilateral Triangles 

The special case of equilateral triangles is particular 
straightforward to interpret because it depends only on one 
scale. Moreover, under the effect of the re-arrangement by 
Eq. (7), an equilateral triangle does not change shape, so 
this case should display self-similar evolution. 

In Fig. (1) the bispectrum of equilateral triangles of 
simulations is plotted against tree-level perturbation the- 
ory, our model Eq. (10) and the SC2001 for reference. The 
rescaled tree-level perturbation theory agrees with simula- 
tions remarkably well at scales k < 0.7h/ Mpc at z — 0. The 
range of scales in agreement also keeps increasing at higher 
redshift, reaching k ~ 3h/Mpc at z = 2. 

In the quasi-nonlinear regime at large scales of fc < 
0.1/i/Mpc at z = 0, Eq. (10) (our model) differs little from 
the tree-level perturbation theory as desired, while the sus- 
picious lowering of SC2001 is clearly seen at low redshifts. 
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Figure 3. B{ki,k2,0) at z=l, (9 = cos ^ /ii2, see Fig. (1) for legend and labels. 



Due to the large cosmic variance, it is not yet possible to 
tell which model is better. 

Our model predicts too much power at scales in the 
strongly nonlinear regime, where the bispectrum of the sim- 
ulations begins to display a power-law behaviour. The break- 
down can be attributed either or both of the two principal 
weaknesses in our consideration: either Eq. (8) is simply too 
rough, or the correlation between window functions at dif- 
ferent values of k is not negligible. 



The major problem of our model for very tilted triangles 
is that if we expand B{k\, k2, 9) with Legendre polynomials 
then it has a very low quadrupole moment (Szapudi 2004). 

The overall performance of Eq. (10), especially at high 
redshift, is comparable with the perturbation theory at one- 
loop level although it is far simpler to implement. In fact, 
as a quick check we realize that the reduced bispectrum Q3 
of our model is very similar to the one-loop results demon- 
strated in Scoccimarro et al. (1998), but emerges in a much 
more straightforward way using our rescaling ansatz. 



3.2 Other Configurations 

The extra complexity of bispectrum over the power spec- 
trum is that it possesses angular dependence as well as 
scale dependence. Comparison of simulations with models 
are drawn in Fig. (2) and (3) for 2 = 0, 1 respectively. 

It is very clear that the our model agrees much better 
with the simulations at higher redshift than the alternatives 
we considered. It is also observed that Eq. (10) is more ac- 
curate when the ratio kijkx < 2; within this range the vari- 
ation with 9 also matches the simulation results rather well. 



4 DISCUSSION 

Following the path pioneered by Hamilton et al. (1991) and 
Peacock & Dodds (1996), we deploy a scale transformation 
argument to construct a well-behaved model of the bispec- 
trum of dark matter fluctuations which is valid in the quasi- 
linear and intermediate nonlinear regimes. On the basis of 
the tests we have been able to perform, the resulting approx- 
imation seems to work exceptionally well, although we are 
somewhat hampered by numerical limitations. A full assess- 
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ment of the precision and reliability of this idea will have to 
wait until simulations of even larger size and higher resolu- 
tion than the VLS are available. Nevertheless, the important 
point may lie even deeper than the pragmatic usefulness of 
this as a simplifying ansatz. The scale transformation of Eq. 
(1) or (2) may contain more physical meaning than has been 
previously thought, and it may not after all be the case that 
higher-order correlation functions in the non-linear regime 
necessarily pose such fierce analytical and numerical chal- 
lenges as has generally been asssumcd. 

The tree-level bispectrum is fully determined by the lin- 
ear power spectrum, while the success of our model seems 
to tell us that the nonlinear bispectrum is governed by the 
nonlinear power spectrum. It therefore seems possible that 
having the power spectrum at hand, one can accurately de- 
rive the bispectrum of a broad range of configurations using 
only the knowledge that it was generated by gravitational 
physics. 

It is well known that a Gaussian random field is fully 
described by its power spectrum, but does a non-Gaussian 
field evolved by gravity from Gaussian initial conditions also 
possess this feature at some level? It will be very interest- 
ing to test this idea by seeing if the trispcctrum can also be 
obtained by rescaling trco-lcvcl perturbation theory. More 
fundamentally, if gravitational interactions do manage to 
act in such a way then why is it that the simple scale trans- 
formation, Eq. (2), manages to capture so much complicated 
physics? 
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